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In the cell, protein complexes form relying on specific interactions 
between their monomers. Excluded volume effects due to molecular 
crowding would lead to correlations between molecules even without 
specific interactions. What is the interplay of these effects in the 
crowded cellular environment? We study dimerization of a model 
homodimer both when the mondimers are free or tethered to each 
other. We consider a structured environment: Two monomers first 
diffuse into a cavity of size L and then fold and bind within the cav- 
ity. The folding and binding are simulated using molecular dynamics 
based on a simplified topology based model. The confinement in the 
cell is described by an effective molecular concentration C ~ L 
A two-state coupled folding and binding behavior is found. We show 
the maximal rate of dimerization occurred at an effective molecular 
concentration C°p ~ ImM which is a relevant cellular concentra- 
tion. In contrast, for tethered chains the rate keeps at a plateau 
when C < C°p but then decreases sharply when C > C°'P. For 
both the free and tethered cases, the simulated variation of the rate 
of dimerization and thermodynamic stability with effective molec- 
ular concentration agrees well with experimental observations. In 
addition, a theoretical argument for the effects of confinement on 
dimerization is also made. 

molecular crow/ding | dimerization | folding and binding | effective molecular 
concentration | Arc homodimer monolayer 
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Many biological functions depend on protein complexes 
or multimeric proteins which must specifically form in 
a crowded cellular environment. There are several types of 
protein complexes. Homodimeric proteins consisting of two 
identical chains or monomers with a symmetrical conforma- 
tion are the most typical 1 . In vitro experiments show that 
the formation of dimeric proteins, termed dimerization, may 
be described as two-state 2 or three-state [3]- Here, the term 
two-state indicates that the folding and binding of monomers 
are directly coupled, while the term three-state signifies that 
binding starts from already folded monomers or that bind- 
ing has a dimeric intermediate. Since dimerization involves 
assembly of two monomers, its rate should depend on the 
monomer concentration. That is, when the separation dis- 
tance of the monomers is large, the monomers should diffuse 
close to each other first and then dimerize. For In vitro ex- 
periments where there is only one kind of molecule involved 
in general, the dimerization occurs easily when the concentra- 
tion is large. In vivo dimerization of specific monomers is more 
complicated than in vitro because cells are rather crowded due 
to the presence of various macromolecules[4][5l[6l[7l|8]. When 
the local concentration of the monomers is low, the monomers 
take a long time to diffuse together, and the diffusion even 
may be kinetically blocked by other molecules. This makes 
dimerization more difficult. Nevertheless, when the local con- 
centration of the monomers is sufficiently high, dimerization 
occurs easily. 

The concentration of total macromolecules in cytoplasm 
is estimated to be 80 ~ 200g/l [4j[5l|6] which is approximately 



ImM (or 100/iM) if the averaged molecular weight m = 
500 X llODa, i.e., 500 amino acids (or 5000 amino acids) in av- 
erage for a macromolecules, is assumed. Obviously, crowding 
must lead to excluded volume effects [71 [51 [9l [TOl [TTl [T2l [T3| 
which can be described using an effective concentration of 
the reacting molecules. Crowding can preferentially destabi- 
lize the balance between reactants and products, and makes 
the association reactions highly favored. It has been sug- 
gested that association constants under crowded conditions 
could be several orders of magnitude larger than those in di- 
lute solutions [5J|5] [71 [S]. At the same time, crowding causes a 
decrease in the diffusion rate of molecules by a factor in the 
range of 3 ~ 10 H [H]. 

The translational diffusion of molecules in the cell is a ki- 
netic process which can be described using Brownian dynam- 
ics |14l 115) . Dimerization, ultimately, involves the intimate 
contact between two specific monomers, a local dynamic pro- 
cess. Previously, the simultaneous folding and binding of a 
number of homodimers has been theoretically studied using 
topology based models (Go-like) by adding a covalent link- 
age between the two monomers of the dimer [161 117] . Such 
studies may be directly related to the in vitro situation. To 
study the in vivo dimerization of homodimeric proteins, both 
interactions between the monomers and those between the 
monomers and other macromolecules must be considered. In- 
cluding crowding effects in such studies may provide some 
useful insights into the formation of various protein complexes 
and protein-protein interactions, and thus enable one to un- 
derstand intracellular protein networks, and to design protein 
complexes that could act as pharmacological inhibitors. 

Here, we study the dimerization of two monomers encap- 
sulated in a cavity with size L which mimics the crowding in 
cell by an effective molecular concentration C ~ L~"^. We 
study both the thermodynamics and kinetics. The diffusion 
of the monomers into the confined space is described by a 
Brownian dynamics. Dimerization depends on the size L, or 
the effective concentration C. There the model predicts a 
maximal rate of dimerization at an optimal confined space 
size = 22 (with unit 3.8A). Such an optimal corre- 
sponds to an effective concentration ~ ImM which is of 
the order of the macromolecular concentration in cells [11[5][H]. 
This suggests a possibility that the rate of dimerization and 
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the concentration of various macromolecules in cells may have 
been optimized by evolution. Based on the changes of the con- 
formational and translation entropies due to the confinement, 
we show that there is a scaling behavior for the heights of free 
energy barriers for binding and the folding transition temper- 
atures with the cavity sizes. 

Results and Discussions 

Molecular Crowding and Molecular Diffusion. Suppose that 
in a cubic box with size L;, = lOOOA, corresponding to a 
small compartment of the cell, there are about 1000 molecules. 
Among these only two specific monomers can form a dimer. 
The effective molecular concentration is Ce — lOOg/l given 
an average molecular weight m — 55KDa (i.e., ~ 500 amino 
acids). That is, 1000 = Ce x L^. All molecules diffuse 
randomly in the large box. Once the distance between two 
monomers reaches a smaller value, L, i.e., the monomers dif- 
fuse into a small confined space (Fig.lA), there they can form 
a dimer by folding and binding. The diffusion time is simu- 
lated using a Brownian dynamics (see the Methods), and a 
monotonic decrease of the time versus the size L is obtained 
(Fig. IB). 

In the Brownian dynamics each monomer is modeled as 
a single particle mimicking the protein chains. Clearly, the 
diffusion of such particle could not fully describe the case of 
the protein chains since the protein chains are soft and may 
change their conformations during the diffusion. However, for 
the sake of simplicity, an effective friction for the particles 
could be used to model the diffusion of the protein chains. 
Previously, a friction coefficient ja ~ 0.05 was used for a single 
amino acid with size a [18]. For the subunit of Arc repressor 
(with 53 amino acids) studied in this work, a corresponding 
friction coefficient is estimated to be 7 ~ 0.19 since the size of 
the Arc monomer is approximately Here, a spherical 

conformation for the monomer is assumed. Note that for a 
protein with 500 amino acids a friction coefficient is 7 ~ 0.4 
[If a friction coefficient 7a — 0.2, which was argued to be a 
factor of 10 larger than the measured value for amino acids in 
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Fig. 1. Moi ecular crowding and confinement. (A) A scliematic particle model for 
the molecular crowding. Two monomers are distributed randomly in a box and dif- 
fuse into a confined space with size L. {B] The time of any two monomers diffusing 
closely with a separating distance < L versus the size L for four different friction 
coefficients. ( C) The Arc dimer (PDB structure larr) confined in a cylindrical cavity. 
{D) The cavity size L versus the effective molecular concentration C ~ 1/{27tL'^) 



water [TQ]. is set, one has 7 ~ 1.6]. Thus a value of 7 ^ 0.1 is 
in the reasonable range to model the kinetics for protein Arc. 
To show the effect of friction coefficient on diffusion rate of the 
monomers, several cases with different 7 values are simulated 
(see Fig. IB). Clearly, the diffusion slows down when the fric- 
tion is large. It is noted that such simplified diffusion of Brow- 
nian particles is approximate but reasonable when the sizes of 
the conformations of the two protein chains can be negligible 
when compared with the inter-chain distance between them. 
Obviously, if the density of the specific monomers is high, the 
diffusion time will be short and the dimerization events will 
be more frequent. In the present work we only put two such 
monomers to model the least case. It is also worth noting 
that if the size of confined space L is roughly 5 ~ 6 times the 
radius of gyration Rg of the monomers, the monomers need to 
diffuse further within the confined space. Here a value 5 ~ 6 
is set since unfolded monomeric chain is extended |llj . 

Model Protein and Confinement. The homodimer studied in 
this work is the Arc repressor of bacteriophage P22 which 
consists of two chains each containing 53 residues and has a 
symmetrical native structure [20] . Clearly, Arc is taken as the 
model protein because it is small and has been studied ex- 
perimentally. In a crowded circumstance, the confined space 
is taken as a cylindrical cavity (Fig. 1(7) and its diameter 2L 
and height h are set equal to each other. The size of the 
cylinder is related to the effective molecular concentration by 
C ~ l/(27rL^) ff every cylinder only contains two molecules. 
Obviously, a big cylinder or weak confinement corresponds to 
a low C and vice versa (Fig. ID). Thus, two Arc monomers 
may fold and bind in such a cavity, modelling the dimeriza- 
tion of a homodimer within a crowded cell. Such behaviors are 
simulated based on the molecular dynamics using the Go-like 
potentials (see the Methods). 

The confinement is modeled by a cylinder which was pre- 
viously used to study the confinement / crowding effect on pro- 
tein folding [5] [15] [13] [H] . Simulations using spherical space 
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Fig. 2. The features of dimerization at within a cavity of L = 20. (A) 
The time evolution of the native contacts: Qj^ for monomeric chain-A in green 
(or Qb for monomeric chain-B in blue), and Qab for the interface in red. {B) 
The time evolution of the separated distance d between the centers of mass of two 
monomers. ( (7) The free energies projected onto versus Qb [D) The free 
energies projected onto versus QaB- 
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provide that using different confinement shapes wiU not qual- 
itatively change the results [221 1101 [23] . Furthermore, a study 
on the molecular crowding effect on folding of globular pro- 
teins suggested that to depict a rather crowded in vtvo envi- 
ronment, the optimal cavity to host protein molecule may be 
cylindrical [11| . 



Two-state Behavior. Some features of the dimerization trajec- 
tories of the Arc dimer confined in a cavity with L = 20 (or 
C ^ 1.2mM) at the related transition temperature are 
shown in Fig. 2. The typical time evolution of the native con- 
tacts of chain-A and chain-B {Qa and Qb) and the interfacial 
native contacts (Qab) and the separating distance d between 
the centers of mass of two chains is shown (Fig.2j4-_B). One 
can see clearly that the folded state (with Qa or Qb ~ 0.9) 
emerges only when the interface is formed (with Qab > 0.85). 
Interestingly, d varies between the folded and unfolded states. 
The free energies of the folding and binding process projected 
onto three different sets of reaction coordinates show the most 
populated states, i.e., folded chains with a well formed inter- 
face (both Qa and Qb ~ 0.9, and Qab ~ 0.85) and the 
unfolded chains without binding (both Qa and Qb ~ 0.5, 
and Qab < 0.1) (Fig.2C-D). Note that the interfacial native 
contacts Nab ~ 143 are almost twice as numerous as the 
intra-chain native contacts Na (or Nb) = 77, thus energeti- 
cally, the states with Qa and Qb ^ 0.5 can still be referred as 
the unfolded states. These indicate that the folding and bind- 
ing occur in a cooperative two-state manner, consistent with 
previous in vitro experimental observation [2l 1241 [25] . and also 
with earlier simulations for linked chains 1161 [TTI . 
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Fig. 3. The features of dimerization. (A) The folding transition temperature 
versus the cavity size L for two free and tethered monomers. {B) The dimeriza- 
tion rate averaged over 100 trajectories at O.SST^ ^ 1.0 versus the cavity size 
L. Curve-A shows the case without the global diffusion for two free monomers, and 
curve-B shows the case with the global diffusion for two free monomers at 7 = 0.1. 
Curve-C shows the case for two tethered monomers. The related concentrations are 
listed in the upper x-axes in both V\g.3A-B. 



Effects of Concentration on Stability. To study the influence 
of the various effective molecular concentrations C on the 
dimerization, the transition temperatures Tf, which charac- 
terizes the thermodynamic stability of the dimer (high value of 
Tf means high stability) are obtained. In Fig. 3^4, it is shown 
that the value of decreases monotonically as L increases 
(or C decreases), implying that small C or large space results 
in low value of Tf or low thermodynamic stability. Exper- 
imentally, both urea and thermal denaturation showed that 
the stability of the Arc dimer is low at low protein concentra- 
tions [21 1241 [^ . Experiments on other dimeric proteins have 
also showed that high concentration improves the thermody- 
namic stability [261 127] . Our results are clearly consistent 
with these experimental observations. The value of Tf at 
C — 22.9mM (or L = 7) increases about 4% with respect to 
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roughly at C = l/iM. Obviously, such a big enhancement in 
the thermodynamic stability is due to the crowding effect or 
confinement which reduces the conformational and transla- 
tional entropies of the unfolded states of the two monomers 
more than it affects the native dimer thus making the unfolded 
states unstable (see an argument in the final part). Note that 
the dimerization cannot occur if the confined space L < 7 (see 
Fig.3A). This relates to a too crowding case for the monomers 
to perform their folding and binding. 



Effects of Concentration on kinetics. The effect of concentra- 
tion on the kinetics of dimerization is also reflected in the 
rate of dimerization by incorporating the diffusion, folding 
and binding processes together (Fig. 35). The rate kf changes 
nonmonotonically as L increases (or C decreases), showing 
an optimal maximum at ~ ImM which is relevant to 
the macromolecular concentration in cells (see curve-A and 
curve-B in Fig. 35). Here the rate kf is in inverse proportion 
to the summation of the time for the two monomers to dif- 
fuse into the conflned space and time for assembly of two 
monomers within the confined space. Note that the assembly 
of two monomers within a confined space may include the lo- 
cal diffusion if the initial distance between the two monomers 
is large. Clearly, here the diffusion of two monomers in the 
cavity is simulated by the motion of two poloymeric chains 
(Fig. 1(7) not of two point particles (Fig.lA). 

In Fig. 35, three cases, namely those for the dimerization 
of two monomers with and without diffusion and for the single 
tethered mutant, are shown. Curve-A shows the case without 
nonlocal diffusion, which describes a situation of high local 
concentration of the monomers. It is found that when L is 
small (or C is high), the dimerization is slow and quite dif- 
ficult since the conformational space for the chains to search 
is limited. As L increases, the dimerization becomes easier 
and faster. However, when L is too large, the conformational 
space becomes very big and the chains now must spend much 
time in finding the folded state, resulting in slow dimerization. 
Thus, there exists an optimal size for the conflnement, or an 
optimal effective concentration C°^. For C < C°^, the rate 
kf monotonically decreases as C decreases (Fig. 35). When C 
is low enough, the rate kf depends linearly on C in agreement 
with the experimental observation [24]. As shown by curve-B, 
a similar change of dimerization rate is also observed when 
the nonlocal diffusion is taken into account. Since the dif- 
fusion time decreases monotonically as the size of conflned 
space L increases (Fig. 15), the decrease of dimerization rate 
becomes slower. However, there still exists an optimal size of 
the conflnement, or an optimal effective concentration having 
about the same value of obtained for the case without 
the nonlocal diffusion. The physical origin for such a behav- 
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ior is basically the same as the case without diffusion, and the 
nonlocal diffusion only increases the total time of the dimer- 
ization when two monomers are further separated. Actually, 
curve-B is related to rather rigorous environment since the lo- 
cal concentration of the monomers is low (only two monomers 
among 1000 molecules are assumed) and the averaged sepa- 
ration distance is large. Clearly, if the local concentration of 
the monomers is not so low or the monomers are co-located, 
the effects of global diffusion are smaller. As a result, a curve 
of dimerization rate should be bounded by the curves - A and 
-B. 

Effect of Confinement for Tethered Mutant. Clearly, for the 
tethered mutant, i.e., when the two chains of the Arc dimer are 
linked together, the thermodynamic stability is higher than for 
the non-tethered case, especially when L > 15 (see Fig. 3^1), 
and the rate of dimerization shows a plateau when C < 
(Fig. 35). Again this agrees with experimental observations 
that tethering two subunits of a dimeric protein significantly 
enhances both the thermal stability of the dimer and the rate 
of dimerization [251128) . The physical reason is that the tether- 
ing reduces significantly the conformational and translational 
entropies of the two tethered chains, resulting in the reduction 
of the search time in the unfolded ensemble and destabiliza- 
tion of the unfolded states. In addition, the two chains of 
Arc dimer would not need to diffuse much to be close to each 
other because they are already linked together. Therefore, 
it takes them less time to complete the folding and binding 
in comparison to the non-tethered case especially for large 
confined spaces. Obviously, more time is needed for diffusion 
as the available space of two monomers grows larger. It is 
worth noting that the tethered case actually is related to a 
rather crowded case of non-tethered monomers, and gives an 
effective concentration Ce = 2.7mM for the non-tethered Arc 
dimers here [25]. This is relevant to the optimized effective 
concentration ~ ImM. 

Free Energy Profiles of Folding and Binding. To character- 
ize the folding and binding of the two chains, we calculate 
the free energy profiles for both processes, respectively. As 
shown in Fig. 4^, for a case of L = 20 (or C ~ 1.2mM), the 
height of free energy barrier AG^ for binding is about 4.7e 
which is much larger than that for folding of the monomeric 
chains, i.e., AG^ ~ 1.7e. This suggests that the binding is 
a dominant rate-limiting step in the dimerization. Interest- 
ingly, it is found that the value of AGj increases when L < 22 
(or G > ImM) and then is saturated to 5.0 when L > 22 
(or G < ImM) (Fig. 45). However, the variation of AG J. for 
the monomeric chain is rather small (Fig.4G). Therefore, the 
crowding effect mainly influences the binding rather than the 
folding of the monomers. The rate limiting step, the binding 
of two monomers, also requires overcoming frustrated polar 
interactions or non-native contacts formed at the interface be- 
tween two monomers in a relatively hydrophobic environment 

To further understand the dimerization of the two chains, 
the free energy profile as a function of Ad, the distance be- 
tween the centers of mass of two chains shifted by subtracting 
the native separation distance, is shown in Fig. 47?. Free energy 
funnels can be clearly seen for three cases. As an example for 
the case of L = 14, a deep well around Ad ~ corresponds to 
a quite stable binding or a localized state of the two chains. 
Note that the effective attraction is short ranged, and is sim- 
ilar to that of the binding between the ligand and receptor. 
It is also seen that the two chains have weak or even no in- 
teraction in a certain range of Ad ~ 12 . However, due to 



the repulsive interaction between the protein chains and the 
cavity wall representing the excluded volume effects of other 
protein, the free energy increases when Ad — Ad* > 16. Thus 
we see the dimerization is cooperatively guided by the binding 
and confinement quite naturally. For the various cavity sizes, 
the ranges with weak interactions and the values of Ad* are 
different, indicating that the slopes of the free energy profiles 
are different. The presence of a free energy funnel allows the 
dimerization to be stabilized by confinement. This is very 
similar to that of a free energy of protein-ligand binding ob- 
tained theoretically and also to the force measured for the 
ligand-receptor association and dissociation 30, 31 . 

The dimerization reaction under confinement conditions 
investigated in this study by the native topology-based model 
|32l 1161 117] focuses on the effect of the confinement on the 
configurational and translational entropy [10]. In the cell, con- 
finement and crowding encapsulate the monomers closely and 
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Fig. 4. The free energy profiles and barriers for folding and binding. (A) The free 
energy for the interface as a function of Qab (solid circles), and the free energy for 
the monomeric chain as a function of (open circles). [B) The height of free- 
energy barrier for the dimer AG^ (marked in Fig. 4^4.) versus L. {C) The height of 
free-energy barrier for a monomeric chain (marked in Fig. 4^4) versus L. {D) 

The binding free energy between two monomers as a function of Arf. 
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Fig. 5. Scaling behavior of folding and binding. (A) The scaling of 
ln(Tj^ - T^^ulkyrpbulk ^j^h tethered monomers. The data are 

taken from Fig.3-A (open circles) and the line represents the theoretical argument 
[Tf - T^^ikyrj^bulk ^ 7^-15/4 b^^^ig^ ^^gjgh^ ^Qt yg^^^j^ ^he 

folding temperature of two free monomers confined in space with L. The data 
are taken from Fig. 45 (solid circles). The line in the main graph is a guide to the 
eye. Inset: The changes of barrier heights of the single monomer scaled 
with the size L. The open circles show the data from Fig. 4(7, and the line shows the 
scaling AAG*. L-l^/'^. 
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facilitate binding. It is possible that the cavity has another 
role besides restricting the available volume of protein mo- 
tions and dynamics. For example, other effects can arise due 
to interactions of the protein with the walls of the cavity or 
due to intra- or inter-monomeric non-native interactions. 



Theoretical Interpretation of Confinement.lt is well known 
that the folded state corresponds to a compact conforma- 
tion, while the ensemble of unfolded states has a huge num- 
ber of extended conformations. Thus, confinement primar- 
ily affects the free energy of the unfolded states through the 
conformational entropy. This effect can be quantified based 
on the theory of polymers with excluded volume [S^. From 
the scaling arguments, the conformational entropy cost reads 
[33 El Ea [10] ASS/fcs oc -iV''/'»(a/L)i'5/4 ^ 
notes the unfolded states, 5"^ the conformational entropy, A'^ 
the number of residues (or beads) with size a of the beads 
in a chain, and L the size of the confined space. The ex- 
ponent 15/4 is more generally equal to 3/(3j/ — 1) where 
V = 3/5 is the Flory exponent. In addition, since at the 
folding transition temperature the free energy differences be- 
tween the unfolded states and native state AG = G" — G" for 
both cases with and without confinement are zero, we have 
a relationship between the folding temperatures and the en- 
tropies as T^t"'*(5'6""= _ s^^'fe) = r/'(S'J - 5^) where the 
superscript L and bulk indicate cases with and without the 
confinement. Thus, we have (T/ - T'^^"^){St''"' - 5^"=) = 

T^KSf: - - {St - 5^'*=)] = T^{AS„ - ASu). In gen- 

eral, we have both the conformational and translational parts 
for AS, i.e., AS = AS" + AS\ 

For the tethered case, two monomers actually become a 
"long" single chain, and their contributions of translational 
entropies to AS are cancelled in a first approximation. Only 
their contributions to conformational entropies remain. Thus 
we have (T/ - rf^"'*)/T/ oc -ASS/fcs oc L'^^/*. Since 
the relative shift {Tf - Tf"-"=) is quite a bit smaUer than 
unity, we have (T/ - T'^'^^'yj.buik ^ ^-15/4^ plotted 
in Fig.5A for simulation data of Fig. 3^, a well agreement 
can be seen. For the case of two free monomers, the pro- 
cess of folding and binding, i.e., a process 1 -|- 2 — > 12, in- 
volves the loss of one chain or monomer, and the transla- 
tional entropies correspondingly cannot be cancelled. For 
the unfolded/unbound states with 2 chains, this contribution 
changes with the confined volume V ~ L'^ as 21ogy, whereas 
it is only logV for the folded/bound state. Thus, we have 
(T/- - T^""=)/r6«ifc ^ _iog V _ AS^/kB. The logarithmic 
term explains why in the Tj' versus L plot (Fig.3A), the curve 
for two monomers does not seem to converge towards a plateau 
while that for the tethered case, where only the algebraic term 
is present, the curve does saturate at large value of L. 

It is clear that the transition state between the unfolded 
and folded states is an ensemble with a non vanishing confor- 
mational entropy due to its smaller spatial extension than the 
ensemble of unfolded states. The transition state ensemble 
is less sensitive to confinement, so its conformational entropy 
is not affected very much by confinement. When the system 
is confined at a given temperature, say at T = j^buik^ ^j^^ 
relative positions of free energies of the folded and transition 
states are not affected, while the unfolded state is destabi- 
lized. As a result, one then needs to increase the temperature 
by an amount Tf — y^"'* to reach the folding temperature. 
The transition state is stabilized by an amount proportional 
to T/- - r)""=. Thus, the barrier for folding AG* should de- 
crease linearly with T^. Such an expectation is consistent 
with our simulation data as shown in the main plot in Fig.5-B 



where a linear behavior is observed. The inset which shows 
the difference between the bulk barrier and that at a given 
confinement is thus an indirect way to check the above linear 
relation between barrier for folding and folding temperature 
shift. 

Conclusion 

A model of confinement effects on dimerization of a typical 
homodimeric protein was studied. It was found that both the 
thermodynamics and kinetics of the dimerization are affected 
significantly by the effective molecular concentration charac- 
terized by the size of cavity. The thermodynamic stability of 
the dimer can be enhanced and the dimerization can be ac- 
celerated as the concentration G increases. An optimal value 
of ~ l.OmM is obtained. This value is of the order of the 
concentration of macromolecules actually found in cells. The 
confinement and binding enhance the folding funnel, stabiliz- 
ing the dimerization of two monomers. 

Materials and Methods 

Methods. 

Molecular difTusion. The diffusion of the molecules (i.e., particles) in a box is 
simulated using a Brownian dynamics as TTl^Vi (t) = Fi (t) — ^Vi (t) + Ti (t) 
jl8| . Here, V, V and TTl^ are the velocity, acceleration and mass of the particles, re- 
spectively. The subscript index i runs from 2 = 1 to Z = 2 for two specific particles 
or two monomers of the dimer, and from i = 3 to i = 1000 for all other particles 
in the box. For the sake of simplicity, all particles are taken to be identical. That is, 
all the sizes are equal to approximately as 53"'^'^'^ CI and the masses are Tn^ = 53/72 
since the Arc monomer has 53 amino acids. Here the size and mass of an amino acid 
are CL and TTL, respectively. Fi is the force arising from the interaction between the 
particles. A hard-core repulsive interaction between the particles and also between the 
monomer and particles is set as = (cTq/?")^^ where the hard-core radius of 

particle is CTq = (53)"'"'''^4.0A, and r is the distance between the particles. And 
an attractive interaction with the 12-10 Lennard-Jones (LJ) potential between the 
two monomers is set as V(r) = 5(o-Q/r)"'"^ — 6((7q /r) T is the white 
and Gaussian random force modeling the solvent collision with the standard variance 
related to temperature by {T{t)r{t')} = 6^kBTS{t — t') where fcfl is the 
Boltzmann constant, T is absolute temperature, t is time, and 6{t — t ) is the 
Dirac delta function. Four values of the friction coefficient from 'y = 0.01 to 0.5 
are used in our simulations (see Fig. IB). The temperature is set asT = 300K. The 
time unit T is accordingly altered following the formula applied for an amino acid 18] 
and other details of the simulation process are the same as for the folding and binding 
(see the following subsections). Based on 100 runs of molecular dynamics simulations 
starting from randomly positions of all the particles and monomers in the box, average 
time for the two monomers to diffuse into the confined space with different sizes L 
is obtained. A periodic boundary condition is used to model the whole cell. 

Topology Based Model of the Homodimer. A Go-like potential is used to 
model the interactions within the Arc homodimer. For each monomeric chain the 
interactions include the virtual bonds angles Vf,o^^_fj^g;g, dihedral an- 

gles yi/ihedra!. non-bonded pairs of the Ca atoms V^a„_,,„„a [for de- 
tails see Ref.!36')]. Here the superscript 5 denotes chain-A or chain-B, respec- 
tively. Note that similar non-bonded interactions are also used for the native and 
nonnative contacts between the inter-chain residues. The native contact is defined 
as occurring when the distance between any a pair of non-hydrogen atoms belong- 
ing to two residues is shorter than 5.5Ain the native conformation of the dimer. 
Thus, the monomeric and interfacial native contacts can be defined. In addition, 
the crowding effect introduces a repulsive potential V^(ri) between the residue 
i and the cylindrical wall when their distance Vi is less than (Tq = 4A. Here 
^'(^0 = E^50[(o■o/2r0^ - 2{ao/2r^f + l]e(ao/2 - n) (for 
details see Ref.[9]). 

Simulations for folding and binding. The simulations were carried out using 
Langevin dynamics and leap-frog algorithm 18'*3T. The native Arc dimer is unfolded 
and equilibrated at high temperature, and then the unfolded conformations are taken 
as starting states for the folding simulations. The energy scale 6=1 and time step 
St = 0.005r are used. Here T = -^Z TYIO?' j € is the time scale with the van der 
Waals radius of the residues a = 5A. All the length is scaled by A = S.SA, i.e., the 
bond length between two Ca atoms. A friction coefficient '^a ^ 0.05 is used. The 
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thermodynamic variables (e.g., the free energy = E(_Q^ — T'logW(Q) 

with E(^Q^ and are energy of the system and the density of conforma- 

tions at Q, respectively) are obtained using the weighted histogram analysis method 
(WHAM) 36 . The free energies for a monomeric chain F{Qa} (or F{^Qb) and 
for the chain-chain binding F{Q Ab) can be calculated. 
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